Hybrid projection methods for large-scale inverse problems with mixed Gaussian priors
نویسندگان
چکیده
Abstract When solving ill-posed inverse problems, a good choice of the prior is critical for computation reasonable solution. A common approach to include Gaussian prior, which defined by mean vector and symmetric positive definite covariance matrix, use iterative projection methods solve corresponding regularized problem. However, main challenge many these that matrix must be known fixed (up constant) before starting solution process. In this paper, we develop hybrid problems with mixed priors where convex combination matrices mixing parameter regularization do not need in advance. Such scenarios may arise when data used generate sample (e.g., assimilation) or different are needed capture qualities The proposed based on Golub–Kahan process, an extension generalized bidiagonalization, distinctive feature both weighting can estimated automatically during Furthermore, training available, various data-driven (including those learned kernels) easily incorporated. Numerical examples from tomographic reconstruction demonstrate potential methods.
منابع مشابه
Efficient Gaussian Sampling for Solving Large-Scale Inverse Problems using MCMC Methods
The resolution of many large-scale inverse problems using MCMC methods requires a step of drawing samples from a high dimensional Gaussian distribution. While direct Gaussian sampling techniques, such as those based on Cholesky factorization, induce an excessive numerical complexity and memory requirement, sequential coordinate sampling methods present a low rate of convergence. Based on the re...
متن کاملIterative Projection Methods for Large–scale Nonlinear Eigenvalue Problems
In this presentation we review iterative projection methods for sparse nonlinear eigenvalue problems which have proven to be very efficient. Here the eigenvalue problem is projected to a subspace V of small dimension which yields approximate eigenpairs. If an error tolerance is not met then the search space V is expanded in an iterative way with the aim that some of the eigenvalues of the reduc...
متن کاملGaussian Markov Random Field Priors for Inverse Problems
In this paper, our focus is on the connections between the methods of (quadratic) regularization for inverse problems and Gaussian Markov random field (GMRF) priors for problems in spatial statistics. We begin with the most standard GMRFs defined on a uniform computational grid, which correspond to the oft-used discrete negative-Laplacian regularization matrix. Next, we present a class of GMRFs...
متن کاملRegularization with Randomized SVD for Large-Scale Discrete Inverse Problems
In this paper we propose an algorithm for solving the large-scale discrete ill-conditioned linear problems arising from the discretization of linear or nonlinear inverse problems. The algorithm combines some existing regularization techniques and regularization parameter choice rules with a randomized singular value decomposition (SVD) so that only much smaller-scale systems are needed to solve...
متن کاملProjection methods for large-scale T-Sylvester equations
The matrix Sylvester equation for congruence, or T-Sylvester equation, has recently attracted considerable attention as a consequence of its close relation to palindromic eigenvalue problems. The theory concerning T-Sylvester equations is rather well understood and there are stable and efficient numerical algorithms which solve these equations for smallto medium-sized matrices. However, develop...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
ژورنال
عنوان ژورنال: Inverse Problems
سال: 2021
ISSN: ['0266-5611', '1361-6420']
DOI: https://doi.org/10.1088/1361-6420/abd29d